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Abstract 


The supercritical CO» Brayton cycle is considered a promising energy conversion system for 
Generation IV reactors for its simple layout, compact structure, and high cycle efficiency. 
Mathematical models of four Brayton cycle layouts are developed in this study for different 
reactors to reduce the cost and increase the thermohydraulic performance of nuclear power 
generation to promote the commercialization of nuclear energy. Parametric analysis, 
multi-objective optimizations, and four decision-making methods are applied to obtain each 
Brayton scheme's optimal thermohydraulic and economic indexes. Results show that for the same 
design thermal power scale of reactors, the higher the core's exit temperature, the better the 
Brayton cycle's thermo-economic performance. Among the four-cycle layouts, the recompression 
cycle (RC) has the best overall performance, followed by the simple recuperation cycle (SR) and 
the intercooling cycle (IC), and the worst is the re-heating cycle (RH). However, RH has the 
lowest total cost of investment (Cit) of $1619.85 million, and IC has the lowest levelized cost of 
energy (LCOE) of 0.012$/(kWh). The nuclear Brayton cycle system's overall performance has 
been improved due to optimization. The performance of the Molten Salt Reactor combined with 
the intercooling cycle (MSR-IC) scheme has the greatest improvement, with the net output power 
(Waet), thermal efficiency y, and exergy efficiency (ye) improved by 8.58%, 8.58%, and 11.21% 
respectively. The performance of the Lead-cooled Fast Reactor combined with the simple 
recuperation cycle scheme was optimized to increase Ciot by 27.78%. In comparison, the internal 
rate of return (IRR) increased by only 7.8%, which is not friendly to investors with limited funds. 
For the nuclear Brayton cycle, the Molten Salt Reactor combined with the recompression cycle 
scheme should receive priority, and the Gas-cooled Fast Reactor combined with the re-heating 
cycle scheme should be considered carefully. 
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Nomenclature 


Symbols 

area, m? 

heat capacity, J/K 

cost, $ 

cash flow, $ 

flow channel diameter, m 
hydraulic diameter, m 
friction coefficient 
specific enthalpy, kJ/kg 
exergy destruction, W 
channel length, m 
mass flow rate, kg/s 


pitch, mm 


"uS s o9 + ya 


pressure, kPa 

pressure ratio 

heat energy, W 

specific entropy, kJ/(kg: K) 
thickness, mm 


temperature, K 


volume flow, m?/s 
work, kJ 


Subscripts 
0 ambient conditions 
cold cold fluid 

hot fluid 

exergy 

work fluid 


net 

state point 
inlet 
minimum 
maximum 
outlet 
total 


revenues 


Acronyms 
GFR 
SFR 
LFR 


Greek letters 


expenses 
compressor 


turbine 


Gas-cooled Fast Reactor 
Sodium-cooled Fast Reactor 
Lead-cooled Fast Reactor 
Molten salt reactor 

Simple recuperation cycle 
Recompression cycle 
Re-heating cycle 
Intercooling cycle 

Specific cost 

Levelized cost of energy 
Internal Rate of Return 
Payback period 

Size parameters 

Area per net output power 
Heat exchanger 

High temperature recuperator 
Low temperature recuperator 
Main compressor 
Decision-making method 
Non-dominated sorting genetic algorithm 
Organic Rankine cycle 
Pre-cooler 

Recompressor 

Turbine 

carbon dioxide 


Supercritical Brayton 


cycle 


surface roughness 
efficiency 

kinematic viscosity 
relative roughness 
density 

heat transfer coefficient 


1. Introduction 


The overuse of fossil energy is one of the major causes of global warming [1]. As one of the 
most efficient and clean energy sources, nuclear energy has been criticized for its high investment 
costs [2] and poor safety [3]. However, scholars have not given up their research on nuclear energy, 
and with the introduction of fourth-generation nuclear reactors, nuclear power generation has 
gradually become a research hotspot [4-6]. Most current power conversion systems in power 
plants use the steam Rankine cycle and gas turbine systems. However, the high reactor outlet 
temperatures and pressures limit the efficiency of nuclear power conversion systems [7]. 
Compared with the steam Rankine cycle [8], the supercritical carbon dioxide (S-COz) power cycle 
has the advantages of higher system thermal efficiency and a more compact power system 
equipment structure [9], which makes it well matched with nuclear power applications. However, 
compared to the advantages of supercritical helium power cycles in ultra-high temperature 
systems (above 1000°C), S-CO» power cycles are more suitable for medium-temperature 
operation (450-600?C) in Generation IV reactors [10]. The advantages of a simple layout and 
compact structure can save the investment cost of nuclear power generation systems, which is 
economically more advantageous than other energy types [11]. Simultaneously, the small footprint 
allows high-power-density nuclear power equipment to be assembled in the factory before being 
transported to the construction site by truck or train [12], facilitating the commercial use of 
nuclear power. 

Currently, research on the supercritical carbon dioxide Brayton cycle (S-CO2BC) has focused 
on system performance optimization [13], working fluid selection [14], techno-economic features, 
and their combinations with other thermal energy systems [15-17]. Ahn et al. [7] conducted a 
comparative analysis of S-CO2BC systems under various cycle layouts and obtained that 
recompressed Brayton cycle systems have a higher thermal system efficiency and emergency 
efficiency than other cycle layouts. Kim et al. [18] analyzed the parameters, such as thermal 
system efficiency and net system output power, as evaluation indices for an S-CO» recompressed 
Brayton cycle system using the pinch-point temperature difference method, optimizing the 
operating parameters of the system to determine the optimal operating conditions of the heat 
exchanger and other equipment. Song et al. [19] investigated the combined cycle system of the 
S-CO2 Brayton cycle and organic Rankine cycle (ORC), optimizing the system's operating 
parameters to improve the thermal performance of the combined cycle system significantly. Bian 
et al. [20] investigated and compared the effects of different control valves on the dynamic 
performance of the SCO2BC system in terms of both thermodynamic and safety performance. 
Wang et al. [21] investigated the dynamic response characteristics under S-CO» Brayton cycle 
failure conditions and proposed corresponding contingency measures. 

All the above studies on the S-CO» cycle system use a single objective function for 
evaluation. Since the single objective function cannot fully reflect the performance characteristics 
of the S-CO» cycle system, many scholars have also established multi-objective functions to 
research key evaluation parameters of the cycle system, such as thermal efficiency, net output 
power, and irreversible energy loss of the system. Wang et al. [22] established multi-objective 
functions, such as system thermal efficiency and net system output power, as evaluation 
parameters in the application of SCO2BC in tower solar power systems to optimize various cycle 
systems, such as the simple recuperation, recompression, and pre-compression cycles. BATTISTI 


et al. [23] used the system thermal efficiency and heat transfer coefficient as objective functions in 
the S-CO2BC cycle to find the optimal performance of the system based on the variation of the 
heat source temperature, the mass flow rate of the working fluid on the heat source side, and the 
maximum temperature of the working fluid, respectively. Nami et al. [24] conducted an 
exergoeconomic analysis of a cogeneration system, a combined heat and power (CHP) system, 
including a gas turbine, SCO2BC, and ORC. The average product unit cost of the system was 
reduced by 0.56$/GJ after optimization, using the sum of investment cost, total endemic damage 
cost, and environmental impact cost as the objective function. Mohammad et al. [25] optimized 
the techno-economics of a simple recuperation Brayton cycle based on the genetic algorithm (GA) 
so that the cycle produces 7196 of the maximum power at the optimum point but at a cost rate of 
only 33% of the maximum power state. Li et al. [26] also considered the ecological performance 
of the system. They used the non-dominated sorting genetic algorithm (NSGA-II) to optimize the 
net output power, thermal efficiency, and ecological performance of the Brayton cycle system 
driven by a hybrid of fossil fuel and solar power. 

In multi-objective optimization, due to the conflict between different objectives, optimizing 
one objective is at the cost of deteriorating other objectives, so obtaining a unique optimal solution 
is difficult. Instead, coordination and compromise are made among them so that the overall 
objective is as optimal as possible [27]. The general optimization algorithms all aim to obtain an 
optimal set of solutions called the Pareto frontier [28]. Many methods have been proposed for 
obtaining compromise solutions from the Pareto frontier. Li et al. [29] used the NSGA-II method 
to optimize the proposed Brayton cycle model. In addition, multi-attribute decision-making 
methods (MADM), such as TOPSIS, LINMAP, and Shannon entropy, were used to select the 
optimal system parameters from the Pareto frontier calculated by the NSGA-II method. Arora et al. 
[30] proposed a thermodynamic model of a simple recuperation Brayton system with 
irreversibility based on finite-element thermodynamic analysis. The proposed model was 
optimized using NSGA-II and a multi-objective evolutionary algorithm based on decomposition 
(MOEAD), and the proposed model was optimized by Shannon entropy, LINMAP, Fuzzy, 
Bellman-Zadeh, TOPSIS, and other MADM methods to find the optimal system parameters from 
the Pareto frontier. Rao et al. [31] proposed a new multi-objective optimization method, the Rao 
algorithm. A combined solar Brayton cycle-power system case study also investigated the 
proposed Rao algorithm. The MADM method was used to rank the Pareto optimal solutions based 
on the average rank. Kumar et al. [30] performed a multi-objective optimization of the Brayton 
cycle system model by NSGA-II, determining the optimal design parameters from the Pareto 
frontier using MADM methods such as Shannon entropy, LINMAP, Fuzzy, and TOPSIS. 

Therefore, most of the research on multi-objective optimization of S-CO2BC based on 
nuclear power generation has been limited to thermodynamic performance. In contrast, some 
research has been conducted on its economics and safety. Also, there are few effective methods for 
determining the unique optimal solution for multi-objective optimization. This study develops a 
mathematical model of the Brayton cycle for four different cycle layouts and four Generation IV 
nuclear reactors. NSGA- II is used to optimize two objectives, 7: and LCOE, simultaneously. The 
results of the four decision methods are evaluated using Taylor diagrams to obtain the unique 
optimal solution. Finally, the comprehensive evaluation results of the multi-index are provided. 
The main contribution of this study is as follows: 
€  Asensitivity analysis of seven key parameters for four different layouts of the Brayton cycle 


was carried out to determine the decision variables and their range of values. The influence of 
the design parameters on the thermodynamics and economics of the system was analyzed. 

€ Taylor diagrams were used to evaluate the Pareto compromise solutions determined by the 
four MADM methods and to determine the unique optimal solution. 

€ The GI+TOPSIS method was used to comprehensively evaluate the optimization results by 
comparing 16 schemes comprising four nuclear reactors and four layouts in terms of safety, 
compactness, thermal, and economic performances. The characteristics and final ranking of 


each scheme are also given. Decision-makers can make scheme choices based on demand. 


2. System layouts and assumptions 


This study uses four different layouts of the Brayton cycle for optimization. One of the most 
basic cycle layouts is the simple recuperation cycle (SR). It consists of a turbine (Turb), a heat 
exchanger (HX), a recuperator (RC), a pre-cooler (PC), a compressor, and a generator, as shown in 
Figure 1. The working fluid is heated in the recuperator (state points 2 to 3) after being 
compressed near the critical point (state points 1 to 2) and then reheated by the cooling medium of 
reactor coolant in the heat exchanger (state points 3 to 4). At this point, the working fluid, which 
has a high enthalpy at the heat exchanger outlet, will expand in the turbine to do work and drive 
the generator to produce electricity (state points 4 to 5). The expanded working fluid is then 
cooled on the low-pressure side of the recuperator (state points 5 to 6) and further cooled by 
exchanging heat with cooling water in the pre-cooler before being compressed again (state points 
6 to 1). 
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Fig. 1. (a) Schematic and (b) T-S diagram of a simple recuperation cycle (SR). 

Fig. 2 shows the recompression cycle (RC), which has one more recompressor and 
recuperator than the simple recuperation cycle. In the recompression cycle, the working fluid is 
split at the outlet of the low-temperature recuperator (LTR) (state point 8), part of the fluid enters 
the pre-cooler (state points 2 to 3), and the rest enters the recompressor (state points 8 to 3b), 
finally converging at the high-pressure inlet side of the high-temperature recuperator (state point 
3b). As the specific heat capacity of the fluid on the low-pressure side of the recuperator is less 
than that on the high-pressure side, increasing the mass flow rate of the fluid on the low-pressure 
side can balance the specific heat capacity on both sides of the recuperator, enhancing heat 


recovery and avoiding the pinch point. 
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Fig. 2 (a) Schematic and (b) T-S diagram of the recompression cycle (RC). 

The other two-cycle layouts are the reheating cycle and the intercooling cycle. As shown in 
Fig. 3 and Fig. 4, the working fluid is heated and expanded twice in the reheating cycle (state 
points 4 to 5 and 6 to 7), which means that more expansion work (Wuri and Wir) can be produced 
for the same compression work input (We), resulting in a higher net output work and cycle 
efficiency. The intercooling cycle regulates the minimum pressure (P1) and intermediate pressure 
(Ps) of the cycle using secondary compression and intercooling, thus reducing the required 
compression work (Wei and W-2) to increase the cycle efficiency. 

Several general assumptions are made to simplify the simulation, as follows: 

(1) The S-CO2BC operates in the steady state. 

(2) The heat losses are neglected during each part of the experiment. 

(3) The cooling water used in the cooler is in the environmental state. 

(4) The resistance loss of the pipeline is ignored. 
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Fig. 3. (a) Schematic and (b) T-S diagram of the reheating cycle (RH). 
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Fig. 4. (a) Schematic and (b) T-S diagram of the intercooling cycle (IC). 


3. Methodology description 


A comprehensive multi-index evaluation and optimization method based on designing power 
conversion systems for nuclear power plants can guide decision-makers. As shown in Table 5, the 
method is based on the hierarchical analysis method, which consists of the scheme hierarchy (four 
nuclear reactors and four-cycle layouts), the index hierarchy (consisting of 11 indices at four 
levels), and the target hierarchy (The best scheme can be evaluated based on the actual needs of 
decision-makers). The data in the index hierarchy is obtained from the model layer. Before the 
final evaluation, it is necessary to optimize each scheme to ensure the accuracy of the evaluation. 
The optimization process mainly uses the NSGA- II algorithm to obtain the Pareto optimal 
frontier. Then, four multi-attribute decision methods are used to obtain the compromise schemes. 
Finally, the optimal results of each scheme were obtained by comparing the four compromise 
schemes through Taylor diagrams. The optimal results are evaluated to obtain the optimal 


scheme. 
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Four types of nuclear reactors Four types of cycle layouts 
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Thermodynamic model 
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Index Hierarchy 
Security Compactness 
Natural prevention capability Area per net output power (4PR) 
Natural mitigation capability Size parameters (SP) 


Thermodynamic Techno-economic 
Net output power (Wnet) Total Cost (Cu) 
Thermal efficiency (7) The Specific Cost ($C) 
Exergy efficiency (72)  Levelized Cost Of Electricity (LCOE) 


Internal Rate of Return (IRR) 


Optimization layer 
Determine the objective functions based on the realistic 
demand. Such as the levelized cost of electricity or the 
net output work 
Optimize the design parameters Single-objective 
optimization Multi-objective optimization 
Acquire the optimal results and analyze the performance 
under the best operating conditions 


y 


Decision-making layer 
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coefficient (Ccoef) and standard deviation (Sstd) to 
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Fig. 5. Description of comprehensive multi-index evaluation and optimization method. 
3.1. Fundamental Model 
This section presents the thermodynamic model, the heat exchanger thermal-hydraulic model, 


and the techno-economic model for the nuclear Brayton cycle system. The models are mainly 
constructed by calling the NIST REFPROP database through the simulation program. 
3.1.1. Thermodynamic model 


A thermodynamic model was constructed based on the first and second principles of 


thermodynamics. The energy calculations for the main components for all cycle layouts are given 


in Table 1. In Table 1, O represents the heat transfer flow rate, W represents output/input power, and h 


represents enthalpy. 
Table 1 


Calculation of energy of main components. 


Components 


Energy 


Heat exchanger Og ono ou — hux im) Se (Tr in — Tr ou ) 
Recuperator Qrecup = yg (hs. m bus) -m,, EE = h,a) 


hy EE (hrr ou ds hrs Jus + hy, in 


Turbine 

We y hry, out — hry in 

Ne out = (he ou is Nein Witte t he in 
Compressor 


Wo = Ne out Nein 


For the recompression cycle, the total input power of the main compressor and recompressor 


can be expressed as: 


Wc =m; (Bussi = yc ou ) SR (1) 


Wre = m, (Rena ur lge out ) e E SR) (2) 


Where: SR is the system circulation diversion ratio. 

In the reheating cycle, as shown in Fig. 3, the reactor cooling fluid is divided equally into two 
parts, which are heated in turn to the working fluid and then converge and enter the reactor, and 
the heat transfer flow rate between the two times is expressed as: 


Qj, -m, (h, -h) = 0.5m,c, (Tin EM G) 
Oi; -m,(h, - h,) 2 0.5m,c, (fs T ai] (4) 
AT, = UL ones ae even -2 (5) 


Where: c, is the specific heat capacity of the reactor cooling fluid, and A T; is the 
confluence temperature difference of the reactor cooling fluid. 
The net output work of the S-CO2BC system: 


n T W rur = Wo (6) 
The thermal efficiency of the S-CO2BC system: 
7], = ¡qe fO x 100% (7) 


Exergy is the part of the energy that can do work, and the exergy destruction is inevitable in 
the actual process. Exergy efficiency is an effective method to evaluate the supercritical Brayton 
cycle. The exergy of each state point can be calculated by: 


E, 2 m| (h, -,) - T, (s; 5.) | (8) 


E cO HS) (9) 


Where m denotes the working fluid mass flow rate, hi and s; are the specific enthalpy and 


entropy at each point, and 0 refers to the environment state. Tr is the temperature of the coolant 
that enters HX to provide heat for the cycle. 
The cycle exergy efficiency is defined as: 


UE =W, |En (10) 


3.1.2. Heat exchangers’ thermal-hydraulic model 

PCHEs were chosen for the heat exchange thermal-hydraulic model because they have been 
adopted in many studies on SCO»-based power cycles due to their compact size and ability to 
withstand high pressures. A model similar to that proposed in Dostal's Thesis [32] is used in this 
work. According to the assumptions established in the Ref. [32], it is sufficient to model a pair of 
hot and cold channels rather than the entire core of the PCHE. PCHEs adopt a semi-circular 
straight channel model [33]. Fig. 6 is the structural diagram of a heat exchange unit. 
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Fig. 6. Structure diagram of PCHE. 

The logarithmic mean temperature difference (LMTD) method was used to calculate the heat 
transfer area, as indicated in Eq.. The total heat transfer coefficient was calculated using the 
contributions of the thermal resistances by conduction and convection (Eq.), and the convective 
heat transfer coefficient was calculated using Eq., where D is the hydraulic diameter. The 


calculation formula is shown in Eq., and the Reynolds number is calculated using Eq.. 


Q = UAM, (1) 
1 1 t 1 
E (12) 
U Arnot A O og 
Nuk 
D (13) 
2 
D- ae (14) 
(17a +a) 
2 
Re = a (15) 
H 


For the hydraulic model, only pressure drops due to friction were considered. The 


well-known Darcy-Weisbach equation (Eq.) was applied in each segment. The total pressure drop 


for each channel is given by Eq. (27). 

In Appendix A, the equations used to calculate the Nusselt number (Nu) and the friction 
factor (f) are presented. It is worth noting that in PCHE, the relative roughness © of the pipe is the 
ratio of the surface roughness € ( € —10? is used in this study) [32] to the pipe diameter d. 

| dd y? 
ASES (16) 
3.1.3. Model validation 

This study intends to establish and verify the model through simulation software based on the 
data of the Advanced Burner Test Reactor Preconceptual Design Report provided by Argonne 
National Laboratory in the United States [33]. This experiment is a performance study of the 
Sodium-cooled Fast Reactor coupled with the S-CO» recompression Brayton cycle (RC); the 
reactor cooling medium is metallic sodium. The initial parameters for cycle simulation are given 
in Table 2 and 

Table 3 for validation. The comparison between the model established in this study and the 
experimental parameters in the references is shown in Fig. 7 and Table 4. The errors of entropy 
and temperature at each point of the cycle are guaranteed to be within 3%, which verifies the 
model's accuracy. Therefore, the model can be used for further parameter analysis and 


optimization. 

Table 2 

Setting parameters in the experiment [33]. 
Term Value 
Turbine inlet temperature Tmax(K) 744.95 
Main compressor inlet temperature Tmin (K) 304.4 
Main compressor inlet pressure Pmin(MPa) 31.25 
Main compressor outlet pressure Pmax(MPa) 7.4 
Turbine isentropic efficiency yr 0.934 
Main compressor isentropic efficiency mc 0.891 
Recompressor isentropic efficiency rc 0.875 
Environment temperature To(K) 303.15 
Environment pressure Po(MPa) 0.1 
Shunt ratio SR 0.71 
Reactor power Qr(MWt) 250 
Sodium mass flow m(Kg/s) 1259 
Sodium inlet temperature Tii (K) 761.15 
Sodium outlet temperature Tz; (K) 606.15 


Table 3 
Geometric parameters of PCHE [33]. 
(mm) d(mm) p(mm) L(m) 
Heat eachanger 2.0 2 2.4 1 
Recuperator 2 1.5 2.3 0.6 


Cooler 1.66 2.0 2.4 0.593 


Table 4 


Comparison between the results calculated using the established model and the data from Ref [33]. 


Parameters Reference Calculated Error (96) 
Working fluid (CO2) mass flow rate ms (kg/s) 1377 1371.4 0.4 
Compressor input power Wc(KW) 53151 56114 5.57 
Turbine output power Wr, (KW) 157530 159570 1.29 
Cooler outlet water temperature Twaterout (C) — 35.8 36.5 1.95 
Thermal efficiency 7;(%) 41.6 41.41 0.46 
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Fig. 7. Validation of the constructed S-CO2BC model. 
3.1.4. Techno-economic model 
Four economic indices and two compactness indices have been adopted in this study to 

holistically assess the techno-economics of the nuclear Brayton cycle holistically. Specifically, 
these include the heat transfer area per net output power (APR) and the turbine characteristic size 
parameter (SP), which represents compactness, and Con SC, LCOE, and IRR, which represent 
economy. 

The heat transfer area per net output power (APR), which is the ratio of heat exchange area to 
net output work, is used as an evaluation criterion for the compactness of the heat exchanger. APR 
is defined as: 


APR = (Ay, + Age HA JW, (17) 


The turbine characteristic size parameter (SP) is an index used to evaluate the ease and 
compactness of turbine manufacturing. The smaller the SP, the more compact the turbine structure. 


It is defined as: 
SP = JV, [Anos (18) 


Where Vr denotes the volume flow rate of the working fluid in the turbine, A hisen denotes 
the isentropic enthalpy drop at the inlet and outlet of the working fluid of the turbine. 

The specific cost (SC) represents the average investment cost of the unit power output of the 
power plant and is defined as: 


ON m 


C,, = (C; F Ls + PREIS Tur, j + x CoN inst (20) 


Where Cio: is the total investment cost, considering the equipment cost (including reactor, 
PCHE, compressor, and turbine) and the auxiliary equipment and installation cost. Cis; is the 
auxiliary installation cost index. 

The cost of PCHE is calculated by referring to the calculation method in Dr Dostal and the 
quotation of Heatric (30 $/kg for stainless steel units) [32], calculating the weight of PCHE to 
estimate its cost. The formula for calculating the weight of the PCHE is given in Eqs. and. 


M pcur 7 Vocug Sn’ P (21) 
zd 

=]——— 22 

p^ 8PI (22) 


Where Vpcue denotes the volume of the PCHE, fm represents the actual volume per cubic 
meter of the heat exchanger. The cost calculation method for other parts is shown in Table 5. 
Notably, among the costs calculated by these formulas, the heat exchangers and reactors cost the 
equivalent of 2003 dollar values, while the remaining components cost the equivalent of 1994 
dollar values [32]. This study converts the calculated cost multiplied by the inflation factor to the 
cost at the end of 2022. Inflation data for the US power sector shows that the inflation rates for 
1994 and 2003 to 2022 are 78.29% and 60.78%, respectively. 

Table 5 


Method of calculating the cost of system components. 


Components Capital investment cost function Data Sources 
Reactor Zr Ci O) C, =2838/KW,, [34] 
1 
Turbine C, =479.34m, f a Je B, )(1+exp(0.0367,-54.4)) — [35] 
2317 


in 


Compressors C, =71.1m,,(1/0.92 — .) PRc- In(PRc) [35] 


The levelized cost Of electricity (LCOE) is the cost of generating electricity over the full life 
cycle of the power plant, defined as: 


PVE 


LCOE = ———————— (23) 
8760u- NY -W 
CF, 
PVE =C,,+ Y ZE (24) 
(19r) 
CF,,, = Wye (OM (1+ er)*) (25) 


The internal rate of return (/RR) is defined as: 


ey I Fra Pos o (26) 
~ (1+IRRÝ 


CF, , = Wye, (8760-u-C,(1- ary ) (27) 


rev, net 


The data required to calculate the economic indices are shown in Table 6. 
Table 6 


Assumptions in the economic calculation [36]. 


Plant data 
Plant installation cost Cins: 1.3 
O&M operations OM ($/kWe) 30 
Escalation rate er (96) 3 
Plant degradation rate dr (96) 1 
Electricity market price Ce ($/kWh) 0.06 
Plant lifetime NY (years) 20 
Plant utilization factor u (%) 85 
Discount rate r (96) 5 


32. Multi-objective optimization model 

Compared to single-objective optimization, multi-objective optimization can optimize two or 
more mutually constrained objectives simultaneously. This study used the non-dominant sorting 
genetic algorithm NSGA- II to optimize two contradictory objectives: 7, (thermodynamic 
performance) and LCOE (economic performance). The turbine inlet temperature Tmax, main 
compressor outlet pressure Pmax, turbine efficiency, compressor efficiency 7c, split ratio SR (only 
for recompression cycle), confluence temperature difference ATe (only for reheating cycle) and 
primary compression pressure Ppr (only for intercooling cycle) were selected as the decision 


variables for optimization. Therefore, the multi-objective optimization model can be expressed as: 


max (7, ) = f (Pond ums Nc, SR, AT ^ pr ) 


min(LCOE) = f, (Paes Tans ire SR, AT 


(28) 
ies mS 


The NSGA- II method is based on a genetic algorithm and the Pareto optimal concept. It is a 
commonly used multi-objective optimization algorithm whose basic idea is to generate different 
sets of decision variables (also called individuals) within a certain range and to use these sets of 
variables to calculate thermodynamic and techno-economic metrics. By evaluating and ranking the 
objective functions, NSGA- II can analyze each optimal solution and the suboptimal solutions 
adjacent to them according to the value of the objective function and plot these solutions on a 
Pareto frontal map to help decision-makers find an equilibrium point to achieve multi-objective 
optimization. The main parameters of the NSGA-II are listed in Table 7. These parameters are 
used to control the execution of the algorithm and influence the quality and computational 
efficiency of the solution. 

Table 7 
Input parameters of NSGA-II. 


Parameters Value 
Population size 100 
Maximum generation 120 
Crossover probability 0.8 


Mutation probability 0.2 


Selection function Binary Tournament 

Objective functions LCOE(min), 7 (max) 

Decision variables Pmax, Tmax, 1c, 11, SR, ^ Tct , Por 
3.3 Decision-making methods (DMM) 

After NSGA-II optimization, only an optimal solution set can be obtained. Further, Pareto 
compromise solutions must be obtained to obtain a unique solution. Previous studies have shown a 
large variation in the Pareto compromise solution obtained by different decision methods, making 
it impossible for decision-makers to judge which method obtains the desired results. Therefore, 
this study uses various decision methods, including Shannon entropy, a modified technique for 
order preference by similarity to ideal solution (TOPSIS), linear programming techniques, 
multidimensional analysis preference (LINMAP), and the analytic hierarchy process (AHP) to 
determine four Pareto compromise solutions. Finally, Taylor diagrams are used to evaluate the 
priority of these compromise solutions to obtain the final Pareto optimal compromise solution. 
3.3.1. DMM1: Shannon entropy model 

The basic idea of the Shannon entropy method is to determine the optimal solution by 
calculating the information entropy value of each solution to measure its stochasticity and 
uncertainty. A smaller information entropy value indicates a more deterministic and reliable 
solution and, conversely, a more uncertain and stochastic solution. 


Normalization of the decision-making matrix: 


E, =——,i=1,2,...,n,j =1,2,...,.m (29) 


Where Ej is the objective value, i and j represent the quantities of the scheme and objective 
function. 
The Shannon entropy is defined as: 
H.- : Y E InE (30) 
j In(n) ij ij 


i=l 


Eq. evaluates the weight of the objective. Eq. calculates the matrix of solutions sorted in 
descending order after decision-making. Shannon Entropy point is searched from the W;, which 
ranked first. 


»-(i-2) 1-8) (61) 
= 
W,=E,-w, (32) 


3.3.2. DMM2 and DMM3: Modified TOPSIS and LINMAP model 

The shortest Euclidean distance between the Pareto optimal solution and the longest distance 
with the nadir point identifies the traditional (TOPSIS) points. The limitation of this method is that 
there may be two TOPSIS decision-making points at different locations. For this reason, relative 
entropy is introduced in this study to correct it. 
Positive Z+ = (z+1 , z*2 , ..., z+n) and negative Z- = (z-1 , z-2 , ..., z-n) ideal solutions are 
explained as below: 


z; - max, ,, {R,} z,=min, , , As] , for the higher the better criteria (33) 


Ab 


z, - min, um {R,} Z; -max,,, (R,) , for the lower the better criteria (34) 


Relative entropy distance is calculated as: 


m zt 1-z* 
D =2 z’ ee) eee (35) 
NE ME ELS 
D 2 z Ea) (36) 


The modified TOPSIS point is selected according to the maximum coefficient T: 


D. 
Le (37) 
D, +D, 


Based on the relative entropy, the modified LINMAP point is defined as: 
L, e min(D; ) (38) 


3.3.3. DMM4: AHP model 

Compared to the previous methods, AHP (Analytic Hierarchy Process) is a subjective 
weighting method. It allows for defining the weight of the objective values in decision-making 
according to the decision-maker's preference. The judgment matrix constructed for the Pareto 
frontier is defined as: 


a. =— (39) 


where 1 and j represent the importance of the corresponding indices, and the importance of 
the thermodynamic, tech-economic, and compactness indices are 3, 2, and 1, respectively, for this 
study. The weight matrix of the AHP method is defined as: 


1 1 
"(Eis | ro eo 
The normalization matrix constructed by the AHP method is the same as that of the Shannon 


entropy model, with specific reference to Eq.. The final decision matrix, Eq., is obtained by 
multiplying the two. The AHP point is searched based on the 4;, which ranked first. 


4 -E, 0, (41) 


3.3.4. Taylor diagram evaluation model 

Taylor diagrams [37] are a graphical technique for comparing performance differences 
between multi-objective optimization solutions. Taylor diagrams allow the performance metrics of 
multiple solutions to be plotted in the same coordinate system to visualize the differences between 
them. The performance metrics include root-mean-square difference (R54), correlation coefficient 


(Ccoef), and standard deviation (S;;4). The specific formulas are expressed as: 


[(6-7)-G-5)] (42) 


R,msa = | 


Cy = E (43) 


1 n ES 1 n E 
Ss, = 246-7) Sue, = = Daan): (44) 
i=l i-l 


3.4 Final evaluation model: G1+TOPSIS 
The optimal result for each solution can be obtained after performing multi-objective 


optimization. In order to comprehensively compare the advantages and disadvantages of each 
scheme, in this study, the Gl [38] (Ordinal Relation Analysis Method) +TOPSIS comprehensive 
evaluation model is established. The four reactors and four-cycle layouts are compared in all 
aspects, including safety, thermodynamics, economics, and compactness. The safety problems of 
nuclear reactors are related to aspects that this study has not discussed deeply. This study refers to 
Bertrand F et al.[4] for the safety assessment of Generation IV nuclear reactors. The safety of 
different nuclear reactors is scored in terms of Natural prevention capability for each concept and 
Natural mitigation capability for each concept. Notably, Ref. [4] does not assess the safety of 
lead-cooled fast reactors, and this study combines the evaluation of the safety of Gen IV reactors 
by Kelly [39] and Tucek et al. [6] to finalize the values of indexes C1 and C2. The concept scores 
for prevention and mitigation are shown in Table 8. Table 9 shows all the indexes used for the 
evaluation. The evaluation process can be referred to Fig. 8. 

Table 8 

Reactor concept scores for prevention and mitigation (x and + means multiplying and adding 
indicators of the same box) [4] 


Prevention indicators /concept GFR SFR LFR MSR 
Loss of heat sink 31 31 31 33 
(Unprotected) loss of flow accident / Loss of 
: 21 21 21 11 
coolant accident 
TOP int t-criticali 
margin 9 prompt-criticality x power T 42 32 T 
increase dynamics 
Passive cooling conduction 11 13 11 33 
Prevention total score 1.9 3.9 3.9 5.6 
Cladding loading in Sub-assembly thermo 
. f 33 22 22 13 
mechanical + chemical 
Primary boundary thermal loading 13 13 13 21 
Containment loading mechanical (static)+ 
f ] 23 33 23 11 
mechanical (dynamic) 
Fission product release inventory x 
12 32 32 22 


confinement/ retention 
Mitigation total score 6.3 7.4 6.8 4.8 


G1 model 


TOPSIS model 
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Fig. 8. G1+TOPSIS comprehensive evaluation process 

Table 9 

Comprehensive evaluation indices. 
Criterion Indexes 


Security (B1) 


Natural prevention capability( C1) 
Natural mitigation capability (C2) 
Wnet (C3) 

m (C4) 

Me (C5) 

Cror(C6) 

SC (C7) 

LCOE (C8) 

IRR (C9) 

PBP (C10) 

APR (C11) 


Thermodynamic (B2) 


Economy (B3) 


Compactness (B4) 


The TOPSIS model is described in section 3.3.2. The Ordinal Relation Analysis Method (G1) 
uses the weighting between each criterion layer. The following briefly describes the steps in the 
Gl method when determining weights. 

Step 1. Determine the ranking relationship of each index according to the importance of each 
evaluation index in the research object. The ranking determined in this study is B1>B2>B3>B4. 
Their relative importance was 1.6, 1.2, and 1.4, respectively. 

Step 2. Give a comparative judgment of the relative importance between adjacent indexes 


and assign a value: 


p -—— (k=n,n—1,:-*,2) (45) 


where, 7 1s the total number of evaluation indexes. 


Step 3. Calculate the index weights w, , and give the Rational assignment w, subjectively. 


(46) 


w3 =rnw, (k=n,n-l,--,2) 


The set of weights derived by the G1 method is (w, ,W3,7,W, j ; 


The final score for each scheme is the sum of the index score for each criterion level 
multiplied by the weights. 


4. Results and discussion 


This section focuses on the effect of seven key parameters on the thermal economy of a 
Brayton cycle system. The feasibility of four decision methods to obtain Pareto compromise 
solutions is analyzed. The differences in the thermo-economic indexes of each scheme before and 
after optimization are compared. Finally, the characteristics of each scheme are given by 
comprehensive evaluation. The main parameters used for the system simulation analysis are given 
in Table 10 and Table 11. The specific parameters of PCHE can be obtained from 
Table 10 


Main parameters of the reactor. 


GFR[40] SFR[41] LFR[42] MSR[42] 
Power rating (MWt) 2400 2400 2400 2400 
Coolant S-CO; Sodium Lead aia ai 
m system pressure 20 0.1 0.1 0.1 
Core inlet temperature (°C) 485.5 371 479 496 
Core outlet temperature (°C) 650 510 573 581 
Core flow rate (Kg/s) 11708 13580 173600 29000 
Power conversion system S-CO2 S-CO» S-CO2 S-CO2 
Table 11 
Conditions of the simulation. 
Input parameters Value 
Compressor inlet temperature Tmin (K) 304.4 
Compressor inlet pressure Pmin(MPa) 31.25 
Compressor outlet pressure Pmax(MPa) 7.4 
Turbine isentropic efficiency yr [43] 0.9 
Compressor isentropic efficiency zc [43] 0.85 
Environment temperature To(K) 303.15 
Environment pressure Po(MPa)[33] 0.1 
Shunt ratio SR 0.71 


4.1 Parametric Analysis 

Firstly, a single-objective parametric analysis was performed to determine the direction of the 
multi-objective optimization. The parameters on the heat source side of this section are analyzed 
using sodium-cooled fast reactors (SFR). Fig. 9 shows the variation of the four indices 
representing thermodynamics and economy with compressor outlet pressure. As Pmax increases, 
the total cost of the system and the exergy efficiency both show an increasing trend. For RC and 
RH, the system's thermal efficiency increases and then reduces as Pmax increases. The levelized 


electricity cost in the system reduces and then increases as Pmax increases. The reason for this 
trend is that as the pressure increases, the enthalpy of the turbine inlet working fluid increases, 
leading to an increase in output power. However, simultaneously, the pressure-bearing 
requirements of the components increase, the system's total cost increases, and the cost increases 
faster than the system's power output, meaning that there is an intermediate Pmax at which y or 
LCOE of the system is optimized. Fig. 10 shows that as Pmax increases, the thermodynamic 
performance of the system increases, and the LCOE decreases. However, the initial investment 
cost (Cio) also increases, meaning that the higher the system temperature within a certain range, 
the better the obtainable thermal and economic performances, but with a higher initial investment. 
Both the cycle efficiency and the cost of the system are influenced by the efficiency of the turbine 
(yr) and the compressor (7c), so this study analyzes their thermo-economic impact. Fig. 11 and Fig. 
12 show the effect of yc and yr on the thermo-economic performance of the system, and they have 
a similar impact on other system indices. Both system y: and ye increase as they increase, and the 
system's thermodynamic performance is enhanced. The system's total cost increases rapidly after 
yr and yc around 0.85. LCOE shows a decreasing trend and then increases, reaching a minimum at 
yr and yc around 0.85. Not the higher the yr and yc, the better it is. Decision-makers should 


consider this carefully. 
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Fig. 9. Effect of main compressor outlet pressure (Pmax) on system thermo-economic performance 
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Fig. 10. Effect of turbine inlet temperature (Tmax) on system thermo-economic performance 
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Fig. 12. Effect of turbine efficiency (yr) on system thermo-economic performance 

This study investigates specific RC, RH, and IC cycle parameters. Among them, the shunt 
ratio (SR) is an important parameter in the recompression cycle, and its influence on the 
thermo-economic of the cycle is shown in Fig. 13. The effect of SR on the thermo-economic of the 
system is monotonic. As the thermodynamic performance of the system decreases with increasing 
SR, the system's total cost decreases. Nevertheless, the LCOE increases with it, meaning that the 
smaller the SR within a certain range, the better, without considering investment costs. The A T; in 
the reheat cycle affects the operating temperature of the system. Therefore, it is necessary to study 
the AT: effect on the thermal properties of the system. As shown in Fig. 14, the thermodynamic 
performance and the system's economy increase with increasing AT. However, the system's total 
cost also increases. A proper primary compression pressure (Py) can bring lower input power to 
the system. The y: and ye junctions of the system decrease Ppr increases, which indicates that an 
increase in Py is not beneficial to the thermodynamic performance of the system. Also, LCOE 
tends to be unfavorable under the influence of Ppr. 
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Fig. 14. Effect of confluence temperature difference (A Te) on system thermo-economic performance 
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Fig. 15. Effect of primary compression pressure (Ppr) on system thermo-economic performance 
4.2 Multi-objective optimization results 

From the analysis in Section 4.1, the cycle parameters Pmax and Tmax affect the system's 
thermodynamics and economy greatly, and the trends of the effects conflict with each other. This 
section conducts a multi-objective optimization study with these parameters as decision variables. 
The constraint conditions of the decision parameters of the cycle with multi-objective 
optimization are shown in Table 12, respectively. Fig. 16 depicts the Pareto optimal frontier of a 
supercritical CO» Brayton cycle system using NSGA-II optimization for four nuclear reactors 
combined with four-cycle layouts. The optimization objectives are the net cycle efficiency (7x) and 
the plant's levelized cost of electricity (LCOE). The LCOE increases moderately and then rapidly 
with the increase of y: because the system-levelized cost of electricity increases faster than the net 
cycle efficiency at higher evaporation pressure and turbine inlet temperature. Similar results have 
been found in previous studies, where thermodynamic and techno-economic performance are 
mutually constrained, with higher thermodynamic performance leading to decreased 
techno-economic performance. Comparing the optimization results of different nuclear reactors, 
the highest y: is for Gas-cooled Fast Reactors, where the highest is for recompression cycles (RC), 
up to 54.67%, due to the higher core exit temperature of the GFR for the same thermal power of 
the reactor, resulting in a higher cycle efficiency. The LCOE of the intercooling cycle is the lowest 
at 0.1207$/KWh, meaning that no single reactor type and cycle layout can achieve the 
techno-economic and thermodynamic performance optimum. The choice can be made according 
to the different needs. 

The ideal point of the Pareto frontier diagram is in the lower right corner outside the 
rectangular frontier, representing the theoretical goal of simultaneously obtaining the maximum 7; 
and the minimum LCOE. To determine the Pareto compromise solution, four decision methods, 
Shannon entropy, modified TOPSIS, LINMAP, and AHP, are applied in this study and compared 
on the Taylor diagram. Using Taylor diagrams, the final Pareto optimal solution is determined 
from these four decision points by measuring the root-mean-square difference, correlation 
coefficient, and standard deviation. Fig. 16 shows the positions of the optimal solutions 
determined by the different decision methods in the Pareto frontier. For the decision points, the 


Shannon entropy point is located at the maximum m point at the top of the frontier, while the 
LINMAP point is in the middle of the frontier. The TOPSIS and AHP points are close to the 


minimum LCOE points, respectively. The Taylor diagram of the decision process is shown in Fig. 


17. 
Table 12 
Bounds of the decision variables. 
Parameters GFR SFR LFR MSR Layouts 
Compressor outlet pressure 
17-30 17-30 17-30 17-30 
Pmax/MPa 
Turbine inlet temperature 803.15~90  663.15-76  726.15-82  734.15-83 All 
Tmax/K 3.15 3.15 6.15 4.15 
Compressor efficiency yc 0.6-0.9 0.6-0.9 0.6~0.9 0.6~0.9 
Turbine efficiency yr 0.65~0.9 0.65~0.9 0.65~0.9 0.65~0.9 
Shunt ratio SR 0.5~0.9 0.5~0.9 0.5~0.9 0.5~0.9 RC 
Confluence temperature 
. 0-10 0-10 0-10 0-10 RH 
difference ATo/K 
Primary compression 
8-12 8-12 8-12 8-12 IC 
pressure P5/MPa 
GFR Pareto Frontier 
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Fig. 16. Pareto optimal frontier and four decision-making points. 
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Fig. 17. Taylor diagram for measuring four decision-making points. 

Fig. 17 shows the scatter distribution of the four decision points in the Taylor diagram. In the 
figure, the correlation coefficient (Cc), root-mean-squared deviation (Rimsa), and standard 
deviation (Ssa) are represented by blue dashed lines, green dashed lines, and black arcs, 
respectively. For example, the Ccoef, Rrmsa, and Ssra of the TOPSIS point of RC in the SFR system 
are 0.9596, 0.0093, and 0.0094, respectively. The values of the ideal point are 1, 0, and 0.0179, 
respectively. The ideal point in the Taylor diagram has maximum Wnet, m1, $e, and IRR, and 
minimum LCOE, Ctot, SC, APR, and SP, which differs from the interpretation of the Pareto 
optimal frontier. These decision points are judged according to the principle that the smaller the 
Rimsa is, the closer the absolute value of Ceveris to 1, and the decision point closer to the ideal point 
is selected. The Shannon point is noticeably always far away from the ideal point, which indicates 
that the single-objective optimal result (maximum m) should not be chosen as the optimized 
working condition. From Fig. 17, the Rimsa of ideal, Shannon entropy, LINMAP, TOPSIS and AHP 
points in the SFR-RC system are 0, 0.0386, 0.0117, 0.0093 and 0.0054, the Cover are 1, -0.9648, 
0.9682, 0.9596 and 0.9529, and the Ssa are 0.0178, 0.0210, 0.0065, 0.0094, and 0.00170, 
respectively. Therefore, the AHP is determined as the final optimal Pareto solution closest to the 
ideal point. Similarly, the Rimsa of the four decision points in the GFR-SR system are 0.01376 and 
0.00417, 0.00512 and 0.0041, and the Cocos are -0.6512, 0.8969, 0.5718 and 0.8245. The Rusa of 
Shannon entropy and TOPSIS points are larger, and the Rimsa of LINMAP and AHP are close, but 
the Cover of LINMAP is closer to 1. Therefore, the LINMAP point is chosen as the optimal point. 


Fig. 18 summarizes the specific values in the Taylor diagram for all scenarios. 


GFR-RC Ccoef Rrmsd Sstd GFR-SR Ccoef Rrmsd Sstd 

Ideal point l 0 0.01815 || Ideal point 1 0 0.00536 
Shannon Entropy -0.46126 0.02470 0.01036 || Shannon Entropy -0.80636 0.01260 0.00787 
LINMAP 0.90497 0.01329 0.00561 || LINMAP 0.87167 0.00352 0.00233 
TOPSIS 0.96124 0.00821 0.01094 || TOPSIS 0.75507 0.00352 0.00381 
AHP 0.95946 0.00529 0.01607 || AHP 0.74271 | 0.00359 0.00402 


GFR-RH Ccoef Rrmsd — Sstd GFR-IC Ccoef Rrmsd — Sstd 

Ideal point 1 0 0.00741 || Ideal point l 0 0.00498 
Shannon Entropy -0.97124 0.02190 0.01463 || Shannon Entropy -0.80800 0.01208 0.00770 
LINMAP 0.98521 0.00420 0.00330 || LINMAP 0.89345 0.00313 0.00226 | 
TOPSIS 0.97763 0.00278 0.00493 || TOPSIS 0.77369 0.00317 0.00358 
AHP 0.95941 0.00212 0.00674 || AHP 0.84020 0.00297 0.00294 | 


SFR-RC Ccoef Rrmsd — Sstd SFR-SR Ccoef Rrmsd — Sstd 

Ideal point 1 0 0.01788 || Ideal point 1 0 0.00620 
Shannon Entropy -0.96479 0.03858 0.02104 || Shannon Entropy -0.65116 0.01376 0.00890 
LINMAP 0.96822 0.01171 0.00649 || LINMAP 0.89694 0.00417 0.00241 
TOPSIS 0.95964 0.00927 0.00937 || TOPSIS 0.57176 0.00512 0.00416 
AHP 0.95292 0.00542 0.01701 || AHP 0.82448 0.00408 0.00303 


SFR-RH Ccoef Rrmsd — Sstd SFR-IC Ccoef Rrmsd 
Ideal point l 0 0.00632 || Ideal point 1 0 
Shannon Entropy 0.38405 0.00756 0.00723 || Shannon Entropy -0.78569 0.01293 
LINMAP 0.25135 0.00613 0.00114 || LINMAP 0.90371 0.00345 
TOPSIS -0.68619 0.01272 0.00752 || TOPSIS 0.67837 0.00419 
AHP 0.90894 0.00339 0.00361 || AHP 0.83549 0.00332 


MSR-RC Ccoef Rrnsd — Sstd | MSR-SR Ccoef Rrnsd — Sstd 

Ideal point 1 0 0.01246 || Ideal point 1 0 0.00549 
Shannon Entropy -0.94002 0.03378 0.02180 || Shannon Entropy -0.77251 0.01287 0.00815 
LINMAP 0.88376 0.00910 0.00403 || LINMAP 0.89387 0.00364 0.00223 
TOPSIS 0.90667 0.00663 0.00726 || TOPSIS 0.77411 0.00354 0.00360 
AHP 0.88784 0.00628 0.00849 || AHP 0.69337 0.00401 0.00444 


MSR-RH Ccoef Rrnsd — Sstd MSR-IC Ccoef Rrnsd — Sstd 

Ideal point 1 0 0.00423 || Ideal point 1 0 0.00538 
Shannon Entropy -0.45376 0.00804 0.00519 || Shannon Entropy -0.77547 0.01250 0.00786 
LINMAP 0.59332 0.00350 0.00169 || LINMAP 0.80172 0.00365 0.00258 
TOPSIS 0.32410 0.00440 0.00321 || TOPSIS 0.79215 0.00332 0.00379 
AHP 0.45403 0.00400 0.00327 || AHP 0.88319 0.00323 0.00274 


LFR-RC Ccoef Rrmsd — Sstd LFR-SR Ccoef Rrmsd 
Ideal point l 0 0.01122 || Ideal point l 0 
Shannon Entropy -0.91707 0.02390 0.01319 || Shannon Entropy -0.79828 0.01309 
LINMAP 0.80257 0.00886 0.00321 || LINMAP 0.88675 0.00372 
TOPSIS 0.92982 0.00476 0.00806 || TOPSIS 0.79049 0.00367 
AHP 0.89692 0.00497 0.01039 || AHP 0.70177 0.00426 


LFR-RH Ccoef Rrmsd — Sstd LFR-IC Ccoef Rrmsd — Sstd 

Ideal point l 0 0.00371 || Ideal point l 0 0.00559 
Shannon Entropy -0.56245 0.00764 0.00491 || Shannon Entropy -0.81545 0.01252 0.00754 
LINMAP 0.89114 0.00257 0.00136 || LINMAP 0.81495 0.00363 0.00293 
TOPSIS 0.83368 0.00245 0.00174 | TOPSIS 0.81696 0.00329 0.00393 
AHP -0.44066 0.00545 0.00268 || AHP -0.24426 0.00606 0.00134 


Fig. 18. Cooef, Rimsa, and Ss of the Taylor diagram. 

Table 13 shows the specific values of the multi-objective optimization results. Notably, 
the shunt ratio (SR) and the primary compression pressure (P,) are concentrated around 50% 
and 8.0 MPa, which are the minimum values of their iteration ranges. The smaller cycle shunt 
ratio improves the narrow point temperature difference of the cycle, increasing the heat 
transfer efficiency of the cycle, which in turn improves the net cycle efficiency and reduces 
the power generation cost. Smaller primary compression pressure can reduce the primary 
compression energy consumption and thus improve the cycle efficiency. The cycle maximum 


temperature is close to the maximum of their iterations, indicating that within a certain range, 
the higher the turbine inlet temperature is, the higher the net cycle efficiency is and the lower 
the LCOE is. The optimal value of the cycle maximum pressure (Pmax) is around 20 MPa or 
30 MPa near the upper limit of the iterative range. A higher pressure does not imply better 
thermodynamic efficiency of the system and the cost of power generation. RC and RH cycles 
can achieve better thermodynamic and techno-economic performance with a lower cycle 
pressure. The yc and yr are in the range of 84% to 89%, which indicates that although the 
efficiency of the cycle is higher with larger yc and yr, the cost of the system also increases. 
An appropriate value can make the thermodynamic and economic performance of the system 
reach the optimal value. Most of the final decision points fall under the AHP, TOPSIS, and 
LINMAP schemes, and the Shannon entropy method does not apply in the optimization 
decision. 

Table 13 

Pareto optimal working conditions and decision solutions. 


Schemes = Pmax/MPa  Tmax/K nelY% n%  SRIV% ATK  P¿/MPa Solutions 


GFR-RC 20.7331 903.1488 8454 8498 50.334  — = AHP 
SFR-RC 18.1709 763.1468 85.19 85.19 50.01 — = AHP 
LFR-RC 19.8486 826.1500 86.14 87.06 50.11  — — TOPSIS 
MSR-RC 20.0013 834.1500 85.22 85.76 50.00 — = AHP 
GFR-SR 29.9826 903.1500 84.49 87.27 — = = TOPSIS 
SFR-SR 29.3426 763.1499 85.44 $88.03 — — = AHP 
LFR-SR 29.9884 826.1500 86.39 88.060 — => == LINMAP 
MSR-SR 29.9921 834.1255 85.1 87.43  — = = TOPSIS 
GFR-RH 20.1769 903.1485 83.78 86.19 — 7.4083 — AHP 
SFR-RH 18.8858 763.1500 87.19 88.75 — 9.7056 — LINMAP 
LFR-RH 20.6056 826.1500 87.07 87.88 — 9.9534 — TOPSIS 
MSR-RH 20.6007 834.1500 86.75 88.65  — 9.998] — LINMAP 
GFR-IC 29.9999 903.1497 84.35 87.56 — = 8.0000 AHP 
SFR-IC 29.9960 763.1494 84.25 87.6 = a 8.0083 AHP 
LFR-IC . 29.8960 826.1494 8444 87.08  — a 8.0060 TOPSIS 
MSR-IC — 30.0000 834.1500 84,33 87.97  — ca 8.0126 AHP 

Table 14 

Results for schemes before optimization. 

Schemes Wet n Ne Crot SC LCOE IRR APR SP 

Unit kW % % 10^6$ $/W $/(KW-h) 96 m/w 

GFR-RC 1135.50 47.40 56.48 1518.84 1738.86 0.014 24 0.513 0.039 
SFR-RC 971.80 40.41 51.81 1574.80 2106.66 0.016 19.3 0.522 0.038 

LFR-RC 1048.80 43.81 54.05 1546.38 1916.76 0.015 21.5 0.467 0.038 

MSR-RC 1052.67 44.21 54.32 1540.87 1902.91 0.015 21.7 0.459 0.039 

GFR-SR 953.05 39.78 45.86 1399.38 1908.80 0.014 24.1 0.223 0.040 
SFR-SR 817.71 34.00 41.23 1428.22 2270.59 0.016 19.9 0.215 0.038 
LFR-SR 880.19 36.76 43.43 1413.58 2087.78 0.015 21.8 0.199 0.039 

MSR-SR 883.33 37.09 43.69 1410.29 2075.53 0.015 22 0.96 0.039 


GFR-RH 702.84 29.34 34.94 1336.55 2472.11 0.017 18 0.459 0.057 


SFR-RH 663.68 | 27.60 34.39 1384.00 2710.96 0.018 16.1 0.335 0.055 


LFR-RH 716.05 29.91 36.22 1370.82 2488.75 0.017 17.8 0.255 0.056 

MSR-RH 718.97 30.19 36.44 1367.98 2473.50 0.017 18 | 0247 0.056 
GFR-IC 957.81 39.98 44.10 1371.33 1861.26 0.013 24.8 0.494 0.040 
SFR-IC — 833.75 34.67 39.52 1395.37 2175.70 0.015 20.9 0.549 0.038 
LFR-IC 890.58 | 37.20 41.69 1383.15 2019.01 0.014 22.7 0.505 0.039 
MSR-IC 893.03 37.50 41.96 1380.29 2009.33 0.014 22.8 0.498 0.039 

Table 15 

Results for schemes after optimization. 

Schemes We n Ne Crot SC LCOE IRR APR SP 

Unit kW % % 10^6$ $/W $(KW-h) 96 m/w 


GFR-RC 1208.81 50.46 62.81 1797.04 1486.62 0.012 28.5 0.263 0.040 
SFR-RC 1012.70 42.11 56.71 1840.04 1816.96 0.014 22.9 0.196 0.039 
LFR-RC 1138.21 47.54 60.59 1905.22 1673.88 0.014 25 0.277 0.039 
MSR-RC 1123.27 47.17 60.56 1841.93 1639.79 0.013 25.6 0.290 0.039 
GFR-SR 990.18 41.33 49.08 1727.39 1744.53 0.013 26.6 0.153 0.035 
SFR-SR 846.09 35.18 44.31 1801.34 2129.00 0.015 21.4 0.156 0.034 
LFR-SR 922.60 38.53 46.85 1806.26 1957.79 0.014 23.5 0.140 0.035 
MSR-SR 916.01 38.47 46.89 1753.89 1914.71 0.014 24.1 0.140 0.035 
GFR-RH 757.27 31.61 37.73 1619.85 2139.06 0.015 21.3 0.279 0.057 
SFR-RH 666.85 27.73 34.59 1734.03 2600.31 0.018 16.9 0.227 0.055 
LFR-RH 718.78 30.02 36.31 1694.69 2357.73 0.016 19 0.254 0.056 
MSR-RH 726.96 30.53 36.73 1714.08 2357.88 0.016 19 0.242 0.056 
GFR-IC 1038.81 43.36 48.81 1713.79 1649.76 0.012 28.3 0.337 0.036 
SFR-IC 897.88 37,33 43.94 1744.44 1942.84 0.014 23.7 0.389 0.034 
LFR-IC 959.08 40.06 46.10 1713.28 1786.38 0.013 25.9 0.356 0.035 

MSR-IC 969.69 40.72 46.66 1740.60 1795.00 0.013 25.8 0.346 0.035 

The specific data for schemes before and after optimization are given in Table 14 and Table 
15. Fig. 19 shows the percentage change of each index before and after optimization. Comparing 
the data before and after optimization shows that the entire system performance improved 
significantly after NSGA-II optimization. The positive change indicates an increase, while the 
negative change represents a decrease. The net output power (Wnet), net efficiency (71), and exergy 
efficiency (ye) all show an increasing trend, which indicates that the thermodynamic performance 
of the system has a great improvement after the optimization of operating parameters. The 
MSR-IC scheme has the most obvious improvement, with the three indices increasing by 8.58%, 
8.58%, and 11.21%, respectively. 

The SC and LCOE show a decreasing trend, and the [RR shows an increasing trend, which 
means that the overall economy of the system is improved. However, the Cio of the system also 
increases. This is because the nt and the LCOE as the optimization objectives do not consider the 
Cto simultaneously, and they are mutually constrained. The LFR-SR scheme has the largest 
increase in Cio, by 27.78%, while its IRR increases by only 7.8%. This optimization improves the 
thermodynamic efficiency of the system and the later revenue of the plant, but the optimization 
results are not friendly to investors with limited funds. For different cycle arrangements, the 
change in SP is not the same, where the compactness of the turbine for /C and SR has respectively 


improved by 1096-1596, while the compactness of RC and RH has decreased, the compactness of 
SFR-RC has decreased by 3.79%, and the rest is not significant. 
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Fig. 19. Percentage change of each index before and after optimization 
4.3 Comprehensive multi-index evaluation results 

After processing the indexes of the optimized scheme using the G1+TOPSIS method, we 
obtained the strengths and weaknesses of each index at each level for each scheme, as shown in 
Fig. 20. Each color line in the diagram represents a combination of reactor and cycle arrangement, 
and the area enclosed by the diagram lines represents the degree of superiority or inferiority of this 
option. From the figure, the area enclosed by RC is the largest, which indicates that RC has the 
best overall performance among the four-cycle layouts, followed by SR and IC, and the worst is 
RH. The C10 index of SR is more prominent, and the C10 index of RC is more depressed, which 
indicates that the heat transfer area required for RC to obtain unit net output power is larger, and 
the system is less compact. Although all the indexes of RH are poor, its C6 index is the most 
prominent among all cycle arrangements, which indicates that RH requires the lowest initial 
investment cost. The C8 index for IC is more prominent, indicating that it has the lowest LCOE. 
Each scheme has different characteristics that decision-makers can choose according to their needs. 
All indexes were assembled to finally obtain the comprehensive evaluation results for each 
scheme, as shown in Fig. 21. The results of the weights calculated in this study using the G1 
method are B1 (0.3972), B2 (0.2482), B3 (0.2069), B4 (0.1478). From the evaluation results, the 
final reactor ranking is MSR>LFR>SFR>GFR, and the ranking for the cycle layouts is 
RC>SR>IC>RH. For the nuclear Brayton cycle, MSR-RC is the best scheme, followed by 
MSR-SR and GFR-RC, and the worst evaluation results are GFR-SR and GFR-RH. The results 
analysis is because safety is considered more important, and MSR has the best safety. The results 


of the cycle layouts ranking are consistent with the previous single-objective analysis. 


C10 


C9 


C7 C6 C7 C6 
|) GFR+SR | GFR+RC GFR+RH GFR+IC |_| SFR+SR SFR+RC |__| SFR+RH SFR+IC 


C10 


C9 


CT C6 CT C6 
M] LFR+SR LFR+RC [| LFR+RH LFR+IC |_| MSR+SR MSR+RC | MSR+RH MSR+IC 


Fig. 20. Comparison of all indexes of each optimal scheme. 
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Fig. 21. Comprehensive evaluation results 


5. Conclusions 


In this study, a mathematical model of the Brayton cycle was developed for four different 
cycle layouts. Parametric analysis was conducted to investigate seven crucial parameters that 
affect the thermo-economic performance of the Brayton cycle. The NSGA- II method was used to 
simultaneously optimize maximizing m and minimizing LCOE. Four decision methods obtained 
the compromise scheme from the Pareto optimal frontiers. Taylor diagrams are used to evaluate 
the four compromise solutions to find the optimal scheme. The changes in the indexes of each 
scheme before and after optimization are analyzed and discussed. Finally, the characteristics of 
each index after optimization for all schemes are given. The main conclusions are summarized as 
follows: 

(1) Higher pressure does not lead to better thermodynamic efficiency of the system and 
power generation costs. RC and RH cycles can achieve better thermodynamic and 
techno-economic performance with a lower cycle pressure. 

(2) The overall performance of the nuclear Brayton cycle system has been improved due to 
optimization. The MSR-IC scheme has the most noticeable improvement, with the net output 
power Wet, thermodynamic efficiency y, and exergy efficiency ye improved by 8.58%, 8.58%, 
and 11.21%, respectively. 

(3) The LFR-SR scheme was optimized to increase Cio by 27.78%, while the internal rate of 
return (RR) increased by only 7.8%. This is not friendly to investors with limited funds. 

(4) Among the four layouts, RC has the best overall performance, followed by SR and IC, and the 
worst is RH. However, RH has the lowest Cot, and IC has the lowest LCOE. 

(5) Considering all indices of the four levels, the final reactor ranking is MSR>LFR>SFR>GFR, 
and the cycle layouts ranking is RC>SR>IC>RH. For the nuclear Brayton cycle, MSR-RC 
should be given priority, while GFR-RH should be considered carefully. 

Supercritical fluid power cycles have great potential for application. In this study, the 
application of the supercritical CO» Brayton cycle in nuclear power plants is considered. The 
constructed mathematical models and optimization decision methods can guide the design of 
nuclear power plants and wider energy fields with potential heat resources. 


Appendix A. Nusselt number and friction factor calculation 
For the semi-circular straight channel PCHE, the Nussle number is calculated by the Gnielinski 
correlation [44]. The Nussle number in Eqs. is calculated Eqs.~. 


Nu = 4.089(Re < 2300) (Al) 


Nutpe-soon 4.089 
5000 — 2300 


Nu = 4.089 + (Re - 2300)(2300,, Re < 5000) (A2) 
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The friction factor (f) used for the Darcy-Weisbach equation (Eq.) depends on the relative 


(Re.. 5000) (A3) 


roughness of the channels (Eq.) and the Reynolds number. The Reynolds number from laminar to 
turbulent flow is calculated by Eqs. ~. 


€ 
se AS 
7 (A5) 
NE 2000,6 < 0.007 us 
^» | 754 exp(0.0065 / 5), 5 > 0.007 
2000, 5 < 0.007 
= (A7) 
t |1160(1/ 8), 6 20.007 
Re, = 2090 (1/8) ^ (A8) 
Re, 2441.196 ^77 (A9) 
The friction factor (f) in Eqs. 1s calculated by Eq.(35) [45]. 
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When Re¡>Re>Rez, the f is obtained from Eq. , where fi is obtained from the formula given by 
Idelchik [46], as shown in Eq. . The Colebrook-White correlation [44] is used to calculate fi, and fn 
is calculated by iteration. When the error between fi and f» is less than 0.01, it can be assumed that 
f 7 fa. Eqs. ~ use the same calculation method. 


f -( - f) exp |- [0.0017 (Re, - Re) + 4 (A12) 


0.032, 6 < 0.007 


zd (A13) 
f oars | LEI as 0.007 
f,=0.11(8 +68/ Re,) ^ (A14) 
F 2 
1 
Jue (A15) 
2log =o) rae 
Siue 3d 
When Re, « Re « Re, 
f, =0.11(6 + 68/ Re,)^ (A16) 
r 2 
1 
Je (A17) 
2log EN PE 
10 Red f , 3 7 
When Re» Re, 
f, 20.11(868/ Re) ^ (A18) 
r 2 
l 
B (A19) 


Í, = 
2log 29 yi 
ài Red f, 3.7 


For more information on the above equations, it is recommended to refer to the references cited in 


this Appendix. 
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